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The centrality dependence of charged-hadron production in relativistic heavy-ion collisions at 
LHC energies is investigated in a nonequilibrium-statistical relativistic diffusion model (RDM) with 
three sources. Theoretical pseudorapidity distributions are compared with preliminary PbPb ALICE 
data at LHC energy of ^/snn = 2.76 TeV for seven centralities, and a previous RDM prediction 
for this energy is confronted with the data. Refined predictions for 5.52 TeV PbPb based on RDM 
parameter extrapolations are presented. 
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I. INTRODUCTION 

First results of charged-hadron production in PbPb 
collisions at LHC energies of ^/snn = 2.76 TeV [212] have 
recently been analyzed in a nonequilibrium-statistical rel- 
ativistic diffusion model (RDM) ^iJ,- So far, the RDM cal- 
culations have been compared to central collision data. 
Now, preliminary ALICE results from the 2010 LHC run 
at different centralities, and for a large range of pseudora- 
pidities —3.7 < rj < 5.1 have become available ^ which 
we use in this work for a comparison with our RDM cal- 
culations, and for a determination of the corresponding 
parameters, which had in only been extrapolated to 
LHC energies. 

In the RDM, the underlying three sources for particle 
production can be traced back to a midrapidity source 
resulting from gluon-gluon collisions, and two forward- 
centered sources arising from valence quark-gluon col- 
lisions. The particle production sources are broadenend 
in time through nonequilibrium-statistical processes such 
as collisions and particle creations, which are described 
based on a Fokker-Planck type transport equation [MZ] ■ 
The additional broadening of the distribution functions 
due to collective expansion then leads to an effective dif- 
fusion coefficient [8]. 

Even in the presence of expansion, drift and diffusion 
coefficients are related through a dissipation-fluctuation 
theorem in the nonlinear version [Si-jT, of the model 
where a nonextensivity parameter q implicitly accounts 
for the effect of long-range interactions that are consid- 
ered to cause the collective behaviour. In the linear 
model, expansion is treated explicitly [8]. Here the re- 
laxation time (or the drift coefficient) and the diffusion 
coefficient are considered as independent parameters. 

An incoherent superposition of the three sources at the 
interaction time - where the integration of the transport 
equation stops - yields good agreement with charged- 
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particle pseudorapidity distributions at RHIC energies. 
It has been shown in [T^l - fT^ within the RDM that at 
RHIC energies of 0.13 TeV (0.2 TeV) the midrapidity 
source generates about 13 % (26 %) of the produced par- 
ticles in a 0-6% central AuAu collision, whereas the bulk 
of the particles is still produced in the two fragmentation 
sources. At SPS, and low RHIC energies of 19.6 GeV the 
effect of the midrapidity source is negligible [l4| . 

At sufficiently high energies (RHIC at 0.13 and 0.2 
TeV), the effect of the Jacobian transformation from ra- 
pidity to pseudorapidity space together with the super- 
position of the sources produces a slight midrapidity dip 
in the charged-hadron pseudorapidity distribution |14| . 
At LHC energies, the effect of the Jacobian for a given 
particle mass tends to be smaller due to the higher aver- 
age transverse momenta |3j, but the dip in the prelimi- 
nary ALICE data is clearly more pronounced [3] than at 
RHIC energies. If confirmed, this would be a strong hint 
towards the significance of the superposition of the un- 
derlying sources, with the fragmentation sources at LHC 
energies moving further apart in pseudorapidity space 
due to the higher beam rapidity. 

The relativistic diffusion model has previously also 
been applied to pp collisions at RHIC and LHC energies 
|15) . Although transport phenomena are not expected to 
be fully developed here due to the small transverse size of 
the system, the RDM yields reasonable results for pseu- 
dorapidity distributions of produced charged hadrons in 
the 77— range where data are available, and provides pre- 
dictions for large values of rj. 

In a recent paper (THj, the TOTEM collaboration has 
indeed reported first experimental pp results at 7 TeV 
and large \r]\. Measured values for dN/dij of charged 
hadrons range from 3.84±0.01(stat)±0.37(syst) at \r]\ = 
5.375 to 2.38± 0.01(stat)±0.21(syst) at |r/| = 6.37. The 
data account for about 95% of the total inelastic cross 
section. A corresponding RDM prediction as shown in 
Fig. 4 of [TF normahzed to 7 TeV CMS NSD pp data at 
midrapidity (which are below inelastic results) turns out 
to be slightly below the TOTEM values at large \r]\. In 
contrast to MC predictions, it has the correct slope. 
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Our model is also suited for asymmetric systems such 
as d+Au and p+Pb, which are more sensitive to the de- 
tails of the distribution functions. At 0.2 TeV we found a 
midrapidity source containing 19 % of the produced par- 
ticles for 0-20% central d+Au collisions JJj . Particle cre- 
ation from a gluon-dominated midrapidity source, inco- 
herently added to the sources related to the valence part 
of the nucleons, had also been proposed by Bialas and 
Czyz [THj . Another three-sources model has recently been 
used to study the influence of the details of hadronization 
and freeze-out on the constituent quark-number scaling 
in the elliptic flow [19| . 

There exist many complementary models which as- 
sume a central source such as the dual parton model 
[5(71 \7T\ , or the quark-gluon string model . Models 
with specific applications to heavy-ion collisions at LHC 
energies have been summarized in [23l [24] . Approaches 
based on gluon saturation such as [25^, ^26] usually do not 
account explicitly for the contributions of the fragmenta- 
tion sources. The RDM provides a simple and transpar- 
ent analytical framework to investigate the interplay of 
central and fragmentation sources in heavy-ion collisions 
at relativistic energies. 

Within the RDM, we investigate in this work the cen- 
trality and energy dependence of the three sources for 
particle production in collisions of symmetric systems at 
RHIC and LHC energies in direct comparison with data, 
and provide a prediction for central collisions at maxi- 
mum LHC energies. 

The energy range considered here for the three-sources 
model covers RHIC energies of ^/sNN — 0.019, 0.062, 
0.13 and 0.2 TeV in AuAu collisions, the presently acces- 
sible LHC energy of 2.76 TeV in PbPb collisions, and the 
maximum LHC energy of 5.52 TeV. 

The ingredients of the model are briefly reconsidered 
in the following section, the calculations of pseudorapid- 
ity distributions of charged hadrons at RHIC and LHC 
energies are discussed in the third section, results and 
the determination of the RDM parameters depending on 
energy and centrality are performed in the fourth section. 
Conclusions are drawn in the fifth section. 



II. RELATIVISTIC DIFFUSION MODEL 

In the Relativistic Diffusion Model, the time evolution 
of the distribution functions is governed by a Fokker- 
Planck equation (FPE) in rapidity space [7] (and refer- 
ences therein) 
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with the rapidity y = 0.5 ■ \n{{E + p)/{E - p)). The 
beam rapidity can also be written as j/bcam = l^ymax ~ 
=F ln(y's]vjv/'7ip). The rapidity diffusion coefficient Dy 
that contains the microscopic physics accounts for the 
broadening of the rapidity distributions. 
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FIG. 1. (Color online) The predicted (dashed g]) RDM 
pseudorapidity distribution function for charged hadrons in 
0-5% central PbPb collisions at LHC energies of 2.76 TeV is 
shown in the upper part of the figure, with the mid-rapidity 
value adjusted to the ALICE data point [2], and RDM pa- 
rameters as in [3] . The solid curve is a x'^-minimization based 
on the three-sources RDM with respect to the preliminary 
ALICE data from 4 that takes the limiting fragmentation 
hypothesis into account, see text. Error bars have been sym- 
metrized using the larger branches. The corresponding RDM 
parameters are given in Table [l] and Fig. 2. In the lower 
part of the figure, calculated pseudorapidity distributions of 
produced charged particles from AuAu collisions (bottom) at 
y/s~NN ~ 0.13 and 0.2 TeV for 0-6% central collisions with 
PHOBOS data |?rH^ are shown for comparison, see 



The drift J{y) determines the shift of the mean rapidi- 
ties towards the central value, and linear and nonlinear 
forms have been discussed [71 [SI [TO]- In the diffusion 
term, a nonlinearity may enter for a nonextensivity co- 
efficient [30] 1 < q < 1.5 that implicitly accounts for 
long-range interactions causing collective expansion. In 
normal Boltzmann statistics one has q = 1, and the dif- 
fusion term is linear ^31 as in Brownian motion. 

The rapidity distribution of produced particles is cal- 
culated from an incoherent superposition of the beam- 
like fragmentation components at larger rapidities arising 
mostly from valence quark-gluon interactions, and a com- 
ponent centered at midrapidity that is essentially due to 
gluon-gluon collisions. All three distributions are broad- 
ened in rapidity space as a consequence of diffusion-like 
processes, and due to collective expansion. 

Here we use the standard linear FPE that corresponds 
to q = I and a linear drift function 



Jiy) = {yeq-y)/Ty 



(2) 



with the rapidity relaxation time Ty, and the equilibrium 
value yeq of the rapidity. For symmetric systems as inves- 
tigated in the present work, y^q = 0, whereas yeq differs 
from in asymmetric systems jl3j . 

The linear formulation of the transport equation corre- 
sponds to the so-called Uhlenbeck-Ornstein [32] process. 
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applied to the relativistic invariant rapidity for the three 
components Rk{y,t) (fc=l,2,3) of the distribution func- 
tion in rapidity space 
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In the Hnear case, a superposition of the distribution 
functions [TT] using the initial conditions Ri^2{y,t — 
0) = Sly ± ymax) with the absolute value of the beam 
rapidities ymax, and R^iy, t ^ 0) ^ 5{y - y^q) yields the 
exact solution. The mean values are derived analytically 
from the moments equations as 

<yiA't) >= J/eg[l-exp(-t/Ty)]=F2/man;exp (-^/Tj^) (4) 

for the sources (1) and (2) with the absolute value of the 
beam rapidity Vmax, and y^q for the local equilibrium 
source which is equal to zero only for symmetric sys- 
tems. For asymmetric systems, the midrapidity source is 
moving |13) , and the superposition of the sources is very 
sensitive to the values of the model parameters. 

Both mean values < 2/1.2 > would attain y^q for t— 00, 
whereas for short times they remain between beam and 
equilibrium values. The variances are 



exp(-2i/Ty)], 



(5) 



and the corresponding FWHM-values are obtained from 
Tk = \/81n2 • (Ik since the partial distribution functions 
are Gaussians in rapidity space (but not in pseudorapid- 
ity space). 

For symmetric systems, the midrapidity source Rz{y, t) 
= Rgg{y,t) has mean value zero. Since the width ap- 
proaches equilibrium twice as fast as the mean value, the 
central source comes close to thermal equilibrium with re- 
spect to the variable rapidity during the interaction time 
Tint- Full equilibrium as determined by the temperature 
would be reached for Tint/Ty ^ 1. The fragmentation 
sources do not reach < 2/1,2 >= during the interaction 
time and hence, remain far from thermal distributions 
in rapidity space, and do not fully equilibrate with the 
central source. 

The charged-particle distribution in rapidity space is 
obtained as incoherent superposition of nonequilibrium 
and central ("equilibrium") solutions of ^ 



dNch{y,t = Tint) 

dy 



+NlMy, nut) + Nl^^Rgg{y, nnt) (6) 

with the interaction time Tint (total integration time of 
the differential equation). 



III. PSEUDORAPIDITY DISTRIBUTIONS 

To calculate pseudorapidity distributions which de- 
pend only on the scattering angle 0, we convert the 
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FIG. 2. (Color online) Left column: Dependence of the 
diffusion-model parameters for heavy systems (central AuAu 
at RHIC energies, central PbPb at LHC energies) on the 
center-of-mass energy ^/snn (in GeV) as determined in this 
work, from top to bottom: (a) Quotient of interaction time 
and relaxation time; (b) total number of produced charged 
hadrons; (c) fraction of charged hadrons in the central source; 
(d) width (FWHM) of the peripheral sources including collec- 
tive expansion; (e) width of the midrapidity source. Results 
refer to charged-hadron pseudorapidity distributions, with ex- 
trapolations to the highest LHC energy of 5.52 TeV, vertical 
line. The diamonds are fit values at RHIC energies of 19.6, 
62.4, 130 and 200 GeV, and at the LHC energy of 2.76 TeV. 
Right column: Corresponding parameter dependencies on the 
number of participants A^part in PbPb collisions at y^sjvjv = 
2.76 TeV for seven centralities. 



results from rapidity to rj space, rj = — ln[tan(0/2)], 
through 



dN dN dy , , , , , dTV 
ciri ay arj ay 



(7) 



using the Jacobian 



J{r], {m)/pT) = cosh(r/)- 

[l + ((m)/pT)'+sinh2(,7)]-i/2. (8) 
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The Jacobian depends on both the particle mass m (or 
the mean particle mass (m) if particle identification is not 
available), and the transverse momentum px through the 
term {iti/pt)'^- However, when comparing with measured 
dN/drj distributions of produced charged hadrons, par- 
ticle identification in rapidity space is not usually avail- 
able, and one has to resort to an approximate expression 
for the Jacobian. In previous works (such as |13j ) we 
had approximated the mean mass by the pion mass (or 
a slightly higher value [33]) since most of the produced 
hadrons are pions, and pt by the mean transverse mo- 
mentum < > of pions at midrapidity, which is about 
0.45 GeV/c at 0.2 TeV center-of mass energy, and 0.5 
GeV/c at 2.76 TeV. 

Here we also consider the effect of more massive pro- 
duced mesons where the effect of the Jacobian is more 
pronounced. At LHC energies of 2.76 TeV in PbPb, pre- 
liminary spectra for identified t:~ , K~ and antipro- 
tons are available [33] as functions of centrality. Since the 
particle to antiparticle ratios are very close to one at this 
energy [33] , we take identical spectra for the correspond- 
ing antiparticles, and parametrize the pt— dependence 
with a suitable functional form |35j . The relative par- 
ticle abundances in central (0-5%) PbPb collisions are 
83% for pions, 13% for kaons, and 4% for protons, with 
the pion fraction increasing to 84% for more peripheral 
(50-60%) collisions. The corresponding < pT > —values 
calculated from the spectra in central (0-5%) collisions 
are 0.51, 0.86 and 1.37 GeV/c for tt,K and p with the 

< Pt > —value for pions decreasing to 0.47 GeV/c for 
peripheral (50-60%) colhsions. 

To obtain the Jacobian J — dy/drj at midrapidity 
y = rj = for the total charged-hadron distribution 
in a given centrality class, we add the resulting mul- 
tiplicity densities for the three particle species consid- 
ered here, and their antiparticles, to obtain Jy=o = 
{dN/dTj\,j=Q)/{dN/dy\y^Q). For 2.76 TeV PbPb central 
collisions, this yields Jy=o — 0.856, which is smaller than 
the result obtained from the analytical expression for the 
Jacobian ^ calculated with (m) — mjr, and the mean 
transverse momentum. 

This indicates that one should consider the full 
Pt— distributions for all charged hadrons when calcu- 
lating the Jacobian. Since the result is, however, ap- 
proximately equivalent to the use of effective values 
for the mass < moff > and the transverse momentum 

< Pt,cS > in ([s]), we choose < TOcff >= m^, and deter- 
mine < PT,eS > such that Jy=o for the total distribution 
of charged hadrons is exactly reproduced. Since the ef- 
fect of the Jacobian transformation is most pronounced 
at midrapidity, this minimizes possible deviations from 
the correct values. For each centrality class, this yields 



< PT.cS >= m,rJy=o/ y 1 - Jy=0- (9) 

For central (0-5%) PbPb collisions at 2.76 TeV, this 
amounts to < pr.cS >= 0.21 GeV/c, and for 50-60% 
centrality < pT,off >= 0.19 GeV/c. With these effective 
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FIG. 3. (Color online) Pseudorapidity distributions for pro- 
duced charged hadrons in 2.76 TeV PbPb collisions as func- 
tions of centrality, from bottom to top: 50-60%, 40-50%, 30- 
40%, 20-30%, 10-20%, 5-10%, 0-5%. Calculated RDM dis- 
tributions (solid curves) have been optimized in x^-fits with 
respect to the preliminary ALICE data from [3], and using 
the limiting fragmentation scaling hypothesis in the region of 
large rapidities where no data are available. The centrality- 
dependent parameter values are as shown in the rhs column 
of Fig. 2. 



mean transverse momenta, the corresponding Jacobians 
for each centrality class are calculated, and the full pseu- 
dorapidity distribution functions for produced charged 
hadrons are obtained in the three-sources model accord- 
ing to ([6|,([8|. The RDM parameters are then optimized 
in a x^— minimization with respect to the available data. 
We show the results in the next section (Figs. 3,4). 

At RHIC energies, we had investigated the dependen- 
cies of the diffusion-model parameters on incident en- 
ergy, mass and centrality in [TS] [T31 [3S]. In particular, 
the centrality dependence seen in the RHIC data is ex- 
actly reproduced in the three-sources model [13l[14j. The 
parameters are shown in Table [T] as functions of the cm. 
energy in central collisions of AuAu, and in central PbPb 
at LHG energies. 

A prediction for central PbPb at 2.76 TeV with ex- 
trapolated RDM parameters had been presented in [3^, 
with the midrapidity point adjusted to the ALIGE data 
2J. It is shown by the dashed curve in Fig. [l] Evi- 
dently, the predicted fragmentation-peak position is too 
far from midrapidity compared to the preliminary data 
^J, and the experimental midrapidity dip is more pro- 
nounced than in the prediction. 

Although these features may be somewhat modified 
once final LHC data with better statistics become avail- 
able, we have performed a x^— optimization of the three- 
sources model solutions with respect to the preliminary 
ALICE data in order to determine the modification of 
the parameters from the prediction. 
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TABLE I. Three-sources RDM-parameters Tint/Ty, Fi, 2, Tgg, and Ugg. Nl"^ is the total number of charged particles, Jigg the 
fraction of particles produced in the central source. Results for <C yi,2 ^ a-re calculated from ybcam and Tint / Ty . Values are shown 
for 0-5% PbPb at LHC energies of 2.76 and 5.52 TeV in the lower two lines, with results at 2.76 TeV from a ^^-minimization 
with respect to the preliminary ALICE data 4 , and using limited fragmentation as constraint. See Fig. |2] and text for the 
evolution of the parameters from RHIC to LHC energies. Corresponding parameters for 0-6% AuAu at RHIC energies are 
given for comparison in the upper four lines based on recent PHOBOS results [29]. Parameters at 5.52 TeV denoted by * are 
extrapolated according to Fig. [2] Experimental midrapidity values (last column) are from PHOBOS '29' for \ri\ < 1, 0-6% at 
RHIC energies and from ALICE 2 for \ri\ < 0.5, 0-5% at 2.76 TeV. 



y'SNN ybeam Tint/Ty < 1/1,2 > ^1,2 Fg 

(TeV) 



dN I 

dri ll-O 



0.019 T3.04 0.97 =f1-16 2.83 1704 

0.062 =f4.20 0.89 Tl-72 3.24 2.05 3003 

0.13 =f4.93 0.89 =f2.02 3.43 2.46 4398 

0.20 =f5.36 0.82 =f2.40 3.48 3.28 5315 

2.76 =f7.99 0.87 T3.34 4.99 6.24 17327 



314±23[2g] 
0.07 463±34[29] 
0.13 579±23I2S] 
0.26 655±49 [29] 
0.56 1601±60 [2] 
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FIG. 4. (Color online) Pseudorapidity distribution of charged 
hadrons in 0-5% central PbPb collisions at LHC energies of 
y/SNN = 2.76 TeV, with the underlying theoretical distribu- 
tions. The result (solid curve) is based on the interplay of 
central (gluon-gluon, dashed) and peripheral (valence quarks 
- gluon, dash-dotted) distribution functions. At LHC ener- 
gies, the midrapidity value is mostly determined by particle 
production from gluon-gluon collisions. 
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FIG. 5. (Color online) Predicted pseudorapidity distribu- 
tion of charged hadrons in 0-5% central PbPb collisions at 
LHC energies of ^/sNN = 5.52 TeV with extrapolated RDM- 
parameters from Table 1 (upper curve). The curve below 
refers to the x^— optimization of the RDM distribution func- 
tions with respect to the preliminary ALICE data |3] in cen- 
tral collisions at 2.76 TeVPbPb, as in Fig. 1. Results for 0.13 
and 0.2 TeV central AuAu are shown for comparison (lower 
curves, as in Fig. 1). 



IV. RESULTS AND RDM-PARAMETERS 

A determination of the RDM parameter for PbPb at 
2.76 TeV through a straightforward x^— minimization 
with respect to the preliminary ALICE data does not 
yield satisfactory results due to the limitation of the 
present data to ?7 < 5.1. In particular, it yields very large 
values for the width Tgg (FWHM) of the central source. 



and correspondingly large pseudorapidity density values 
at ry = yboam which violate the limiting fragmentation 
scaling hypothesis [23] . In order to obtain realistic RDM 
parameter values, one definitly needs data in the frag- 
mentation region. These will, however, not be available 
in the next years. 

As a remedy, we use the effect of limiting fragmentation 
scaling [37] as a constraint. It has been observed to hold 
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with considerable accuracy in pp and AA collisions [251 
I29j : At sufficiently high energy, particle production in the 
fragmentation region becomes essentially independent of 
the collision energy. The limiting fragmentation concept 
refers to particle production as function of rapidity y, but 
it also holds in pseudorapidity space since for \r]\ 3> 1 we 
have y ~ In {pt I'^t) +1] rj. (In the large-?] region, the 
Jacobian is very close to unity.) 

To supplement the missing LHC data in the pseudo- 
rapidty region close to ybeam, we consider analogous cen- 
trality classes in the 0.2 TeV AuAu results at RHIC, 
select a small segment of RHIC data (up to five of the 
outmost datapoints of the PHOBOS datasets), and shift 
them by Ay = y^HC _ yRfflC ^ 7.99 _ 5.36 = 2.63. Since 
the centrality classes of the measurements at RHIC by 
PHOBOS (221 and at LHC by ALICE do not match 
exactly for more peripheral collisions, a slight correction 
of the shift is performed, except for the two most central 
classes. 

The RHIC data points in the fragmentation re- 
gion are then taken as an additional constraint of our 
X^— minimization of the RDM parameters with respect 
to the preliminary ALICE data at 2.76 TeV. It turns out 
that in contrast to an unconstrained x^— minimization, 
this yields physically reasonable distribution functions 
and RDM parameters, which can be used as basis for 
an extrapolation to the highest LHC energy of 5.52 TeV. 
The midrapidity value is mainly determined by the cen- 
tral source at LHC energies, and is therefore not much 
affected by the use of limiting fragmentation at LHC en- 
ergies. Small scaling violations which may occur in the 
forward region would essentially be disconnected from 
the midrapidity source, although they can slightly mod- 
ify its weight and width. 

The result of the constrained optimization of the 
RDM-parameters with respect to the preliminary PbPb 
data at 2.76 TeV is shown by the solid curve in Fig. [l] for 
central collisions, with parameters given in Table |lj The 
integration is stopped at the optimum values of Tint/Ty, 
^ 1,2, gg, and Ugg] the explicit value of Tint is not needed. 
The fraction of particles produced in the central source is 
Ugg. In view of the preliminary character of the ALICE 
data, we do not list numerical results here. 

The dependencies of the parameters on the center-of- 
mass energy for central collisions is shown in the Ihs col- 
umn in Fig. [2] and their centrality dependence for 2.76 
TeV PbPb in the rhs column of the same figure. Solid 
curves are fit functions. 

The value of the time parameter Tint/Ty shown in the 
upper frame of Fig. [2] is decisive for the position of the 
fragmentation peak in 77— space. It was found to decrease 
for increasing log ^/snn from AGS to the highest RHIC 
energies with a functional dependence discussed in [3], 
and hence, the extrapolation to LHC energies was based 
on a continued fall, resulting in Tint/Ty ~ 0.67 at 2.76 
TeV. 

The comparison with the preliminary data, however, 
reveals that it actually levels off at LHC energies to a 



value of Tint/Ty ~ 0.87 at 2.76 TeV. This indicates that 
in the large energy gap between the highest RHIC and 
the current LHC energy, the rapidity relaxation time Ty 
decreases faster than the interaction time Tint- 

The total number of produced charged hadrons Ntot, 
and the fraction Ugg produced in the central source are 
also displayed in Fig. [2j The latter reaches about 0.56 
at 2.76 TeV, whereas at RHIC energies, it remains below 
0.3, such that the multiplicity density at midrapidity has 
a substantial contribution from the overlapping fragmen- 
tation sources 3J . To decide whether a saturation of ngg 
is attained at LHC energies, one needs the results from 
charged-hadron production at the LHC design energy of 
5.52 TeV, to be expected in 2014/15. 

The partial widths as functions of energy are displayed 
in the lower two frames of Fig. [2] for both fragmenta- 
tion and midrapidity sources. These widths are effective 
values, because they include the effect of collective ex- 
pansion in addition to the statistical widths that can be 
calculated from a dissipation-fluctuation theorem [6]. 

Should future LHC experiments at large pseudora- 
pidity reveal small limiting fragmentation scaling viola- 
tions, this would lead to a modification of the widths and 
weights of the fragmentation sources, and indirectly also 
to small modifications of the central source. 

The values at RHIC energies are resulting from cor- 
responding x^-minimizations with respect to PHOBOS 
data [55]. Note that RDM parameters at 0.13 and 
0.2 TeV (Table [l are slightly different from the val- 
ues in Table 1 of [3j which were based on previous data 
[27| l28] . The data extend in pseudorapidity up to and 
even sligthly beyond yboami with a sizeable pseudorapid- 
ity density measured a.t rj — ybcam- 

In the rhs column of Fig. [2] the centrality depen- 
dence of the RDM parameters as obtained from the 
X^— minimization including limiting fragmentation is 
shown for 2.76 TeV PbPb, and seven centrality classes as 
indicated by the diamonds. The dependences are shown 
as functions of the average numbers of participants for 
each centrality class, from — 5% central to 50 — 60% 
peripheral collisions (see also Fig. [3]). For all five RDM 
parameters, we find a nearly linear dependence on the 
average number of participants in each class. 

There is only a slight decrease of the time parame- 
ter Tint/Ty and the particle fraction in the midrapidity 
source ngg from central to peripheral collisions. The 
strongest centrality dependence is found for the total 
number of produced charged hadrons, which falls from 
about 17,300 in a — 5% central collision to below 1,800 
in a 50 — 60% peripheral collision. The relative strength, 
and the width of the central source decrease only slightly 
towards peripheral collisions. The width of the fragmen- 
tation sources is almost independent of centrality, and 
slightly smaller than the width of the central source. 

The centrality dependence of the pseudorapidity distri- 
butions of produced charged hadrons in 2.76 TeV PbPb 
is displayed in Fig. [3] with the dependence of the RDM 
parameters on the number of participants as displayed in 
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the rhs column of Fig.[2]for seven centralities as indicated 
in the caption. 

At LHC energies, and for aU centrahties, the overall 
scenario changes in favor of particle production from the 
midrapidity source, as can be seen for central collisions 
in Fig. |4j with the Ugg values for the fraction of particles 
in the midrapidity source given in Table |Tj At 2.76 TeV, 
the bulk of the midrapidity density is generated in the 
central source, there is a relatively small overlap of the 
fragmentation sources at midrapidity. 

The sizeable dip at midrapidity that is seen in the pre- 
liminary ALICE data at 2.76 TeV |5 is likely due to the 
Jacobian plus the interplay of the three sources, with 
the fragmentation sources moving much further apart as 
compared to the highest RHIC energy of 0.2 TeV. The 
Jacobian transformation has an effect only on the midra- 
pidity source, which is flattened and has a slight dip. 
In the three-sources approach, the midrapidity minimum 
seen in the data is then easily achieved. 

There exist detailed microscopic calculations of frag- 
mentation sources from gq q and qg ^ q diagrams by 
Szczurek et al. [381 139 1 for pion production in proton- 
proton, and nucleus- nucleus collisions at SPS and RHIC 
energies. These processes are also responsible for the ob- 
served differences |40| in the production of positively and 
negatively charged hadrons, in particular, pions. Extend- 
ing these calculations to LHC energies would lead to a 
microscopic foundation of our three-sources approach at 
energies beyond RHIC. 

The existence and relevance of the fragmentation 
sources is corroborated by results of a partonic model 
for net-baryon distributions that we had presented and 
discussed in |41H43j . For baryons minus antibaryons (and 
also for net protons) , the midrapidity gluon-gluon source 
cancels out, and only the fragmentation sources remain. 
At high SPS and RHIC energies (and accordingly, in 
the predictions for LHC), these give rise to two pro- 
nounced fragmentation peaks which are clearly seen in 
the data. Note, however, that the fragmentation peaks 
in net baryons occur at larger rapidity values as compared 
to charged hadrons. As an example, the net-baryon peak 
in 2.76 TeV PbPb ist at ?/pcak = 5.7 [M], whereas the 
charged-hadron fragmentation peak is at = 3.3. 

With extrapolations of the time parameter in the rel- 
ativistic diffusion model, the numbers of charged parti- 
cles in the sources, and the partial widths Ti^2,^gg from 
Fig. [2] and Table |l] a prediction for central PbPb at 5.52 
TeV is shown in Fig. [5] Here the midrapidity value has 
been extrapolated with log ^/s~nn, and the value of Tgg 
is determined accordingly. The central collision results 



for AuAu (0.13 and 0.2 TeV) and PbPb at 2.76 TeV are 
displayed for comparison. 

V. CONCLUSION 

We have analyzed recent preliminary ALICE results 
for PbPb collisions at LHC energy of y^s^vw = 2.76 TeV. 
Charged-hadron pseudorapidity distributions have been 
calculated analytically in the non-equilibrium statistical 
relativistic diffusion model RDM. For seven different cen- 
tralities, the underlying RDM parameters have been de- 
termined in a x^-optimization of the analytical model 
solutions with respect to the preliminary data, and using 
limiting fragmentation as an additional constraint. 

A comparison with a previous prediction [3, that 
was based on an extrapolation of the parameters with 
log -y/sjvw reveals that the rapidity relaxation time Ty de- 
creases substantially in the energy region between RHIC 
and LHC energies, leading to a larger time parameter 
Tint /Ty and hence, to a fragmentation-peak position that 
is closer to midrapidity than expected from the earlier 
extrapolation of the time parameter. 

Based on the RDM fit to the data in the three-sources 
model, the midrapidity source that is associated with 
gluon-gluon collisions accounts for about 56% of the to- 
tal charged-particle multiplicity measured by ALICE in 
central PbPb colhsions at 2.76 TeV. 

The fragmentation sources that correspond to parti- 
cles that are mainly generated from valence quark - 
gluon interactions are centered at pseudorapidity values 
< r/1.2 >— < 2/1,2 >— t3.3. The total particle content in 
these sources amounts to about 44% of the total charged 
hadron production, but contributes only marginally to 
the midrapidity yield. It is, however, decisive for the 
more pronounced midrapidity dip at LHC energies where 
the fragmentation sources move much further apart than 
at RHIC energies. 

With the results for PbPb at 2.76 TeV LHC energy 
and previous RDM results for AuAu collisions at RHIC 
energies, we have extrapolated the three-sources model 
parameters to the LHC design energy of 5.52 TeV for 
PbPb, and calculated the corresponding charged-hadron 
pseudorapidity distribution. Small corrections of the 
extrapolated values for the diffusion-model parameters 
may be required once the final measured distributions 
become available at both LHC energies. 

This work has been supported by the ExtreMe Matter 
Institute EMMI. 
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